Application of the Lifshitz theory to poor conductors 
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The Lifshitz formula for the dispersive forces is generalized to the materials, which cannot be 
described with the local dielectric response. Principal nonlocality of poor conductors is related with 
the finite screening length of the penetrating field and the coUisional relaxation; at low temperatures 
the role of collisions plays the Landau damping. The spatial dispersion makes the theory self 
consistent. Our predictions are compared with the recent experiment. It is demonstrated that at 
low temperatures the Casimir-Lifshitz entropy disappears as T in the case of degenerate plasma and 
as for the nondegenerate one. 
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The Casimir-Lifshitz force is a dispersion interac- 
tion of electromagnetic origin acting between neutral 
bodies without permanent polarizations. The original 
Casimir formula [H for the force between ideal metals 
was extended to real materials by Lifshitz and cowork- 
ers [1, H, [3l- Recently there was considerable progress 
in the experimental verification of the theory (see review 
and references therein) and various applications to 
nanosciences were discussed [5]. The Lifshitz theory is 
extensively accepted as a common tool to deal with dis- 
persive forces in physics, biology, chemistry, and technol- 
ogy. 

The zero temperature contribution to the force origi- 
nating from quantum fluctuations of the electromagnetic 
field is well understood. On the contrary, the classical 
or high temperature part of the thermal contribution (no 
dependence on K) is the source of constant controver- 
sies. There is no continuous transition for the forces be- 
tween ideal metals and real metals [6] . It is related with 
the transparency of real metals for s-polarized (trans- 
verse electric) low frequency field. It was found that for 
real metals there is a thermodynamic problem the 
Casimir-Lifshitz entropy is not going to zero at T — > 0. 
This problem is still in debate (for the latest publica- 
tions see Q) but a new controversy for poor conductors 
emerged [!|]. In this case the reflection coefficient for 
p-polarization (transverse magnetic) is discontinuous in 
the transition from zero to arbitrarily small conductivity. 
This discontinuity again breaks the Nernst heat theorem. 
Obvious contradiction to the common sense shows that 
some important physics is missed. 

In this paper we demonstrate that account for the spa- 
tial dispersion of the materials resolves the problems and 
make the theory self consistent. We formulate the con- 
dition at which the Lifshitz formula can be extended to 
the description of the forces between nonlocal materials. 
Special attention is paid to the case of poor conductors 
because the nonlocal effects are more important for them 
than for dielectrics or good metals. 

For dielectrics the nonlocality can be neglected due to 
the absence of free charges except maybe the range near 
polariton resonances. For metals the local approximation 
is good because of the very short Thomas-Fermi screen- 



ing length. Spatial dispersion for metals is important at 
low temperatures when the mean free path for electrons 
becomes larger than the fleld penetration depth and the 
anomalous skin effect plays role [13, El, ■ 

Recently it was demonstrated by Pitaevskii that 
the classical part of the Casimir-Polder force between a 
bad conductor and an atom is essentially nonlocal due to 
the finite screening length of free charges (Debye length). 
In Ref. [l^ the force was found in the large distance 
limit. The formula for the force interpolates between 
good metals and dielectrics. It depends on the density of 
free carriers via the Debye length Id. 

The most general theory allowing the spatial dispersion 
was developed by Barash and Ginzburg The force 
was found to have an additional term in comparison with 
the Lifshitz formula related with the nonlocal material 
response. On this bases it was concluded in Ref. [lB| that 
the Lifshitz formula cannot be used in the nonlocal case. 
However, this conclusion does not follow from The 



Lifshitz formula can be applied at least for plasma-like 
media. It was demonstrated that if the long range 
interaction is taken into account in the plasma dielectric 
functions to the first order in the coupling constant e^, 
then the Lifshitz formula holds true. 

Below it is shown that the dielectric functions of the 
plasmas in metals and semiconductors calculated in the 
random phase approximation (RPA) can be used with the 
Lifshitz formula and give sufficient description of these 
materials for evaluation of the Casimir-Lifshitz force. 

To account for the nonlocal response of metals it is pos- 
sible [nl, im to use the RPA (electrons are independent 
but resp ond not to the external field but to the screened 
one [TjI)- This simple choice is good in the weak-coupling 
regime (rarefied plasma), but the plasmas in real metals 
are actually strongly coupled. However, significant devi- 
ations from RPA appear at large wavenumbers k ^ kp 
[l8| . where kp is the Fermi wavenumber. For the bodies 
separated by the distance a the important wavenumbers 
are k ^ 1/a ^ kp. For this reason the use of the RPA 
dielectric functions is justified. 

We can use the same approach to construct the nonlo- 
cal response of nondegenerate semiconductors. The dif- 
ference with metals is that now electrons (holes) have to 
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be considered as nondegenerate plasma and instead of 
the Fermi-Dirac distribution the Bohzmann distribution 
can be used. Technically the problem is equivalent to the 
weak-coupled nondegenerate plasma, for which one can 
use the textbook result ,19]. 

Generalizing the situation we can consider two classes 
of plasmas: degenerate and nondegenerate. The first 
class (I) describes materials for which density of free 
charged particles stays finite at T = 0. Metals, semimet- 
als, degenerate semiconductors etc. belong to this class. 
Materials with the energy gap, for which density of free 
charges disappears with T belong to the second class (II). 
Representatives of this class are nondegenerate semicon- 
ductors, ionic conductors, many disordered materials etc. 

In the nonlocal case the dielectric function becomes 
a tensor, which has two independent components: 
transversal St and longitudinal si with respect to the wave 
vector k. For very small relaxation frequency 7 — > (col- 
lisionless plasma) in RPA these components are 



et.i{uj,k) = Sq{uj) 



uj{uj + 27) 



/m(^)- 



(1) 



Here the first term £o('-^) is introduced to account for the 
interband transitions since these processes are beyond 
the (quasi)free electron model. The variable x is 
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where m* is the effective mass of the charge carrier (elec- 
tron or hole) and Ui? is the Fermi velocity. The plasma 
frequency ujp in ^ and the Debye wavenumber ko (see 
below) are defined as 



ko = 



12^ for /, 



for //, 



(3) 



where N is the density of carriers. Note that £t and 
El are of the first order in and can be used with the 
Lifshitz formula. The functions ft,i{x) are different for 
each class and can be found in the textbooks (see, for 
example, [ll]). For the class I they are 
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For the class II these functions arc 

ft{x) = -Fix), fi{x) ^ -2x^1 + F{x)), 
where F{x) is defined as 



(4) 



(5) 



The dielectric functions ([T]) can be generalized to the 
case of finite relaxation frequency 7. For the transverse 
component et{uj,k) this is an easy task since the colli- 
sion integral in the relaxation time approximation gives 
the same result as ([T]) but with the finite value of the 
relaxation frequency 7. 

Different procedure has to be used to account for the 
finite relaxation time in ei(uj,k). The longitudinal field 
influences the charge distribution. In this case the relax- 
ation of the perturbed distribution toward "equilibrium" 
will be to the local state of charge imbalance and not 
to the uniform distribution [23] ■ Mermin [l^] used this 
idea to generalize Si{k, oj) calculated in RPA to the finite 
relaxation time: 



ei{uj,k) = eo(^) + 



{l + t^)[ei{uj + ij,k)-eo{cu)] 

-1 1 • 7 £i(t^+i7,fc)-£o(") ' 
^ ^uj e,(0,fc)-eo(0) 



(7) 



where we denote the dielectric functions with the finite 
relaxation frequency as et,i. 

The dielectric functions £i^t{k, w) are defined for an in- 
finite medium. The spatial dispersion close to the body 
surface needs special attention. Strictly speaking one can 
define the nonlocal dielectric function only for infinite 
body, but in special cases of specular and diffuse reflec- 
tion of electrons on the surface of the body it is also 
possible to do. We will consider here the case of spec- 
ular reflection. An electron reflected from the interface 
with vacuum cannot be distinguished from that coming 
from a fictitious medium on the vacuum side. In this 
way the specular condition continues the medium with 
the interface to the infinite medium. 

For actual evaluation of the force in the case of ma- 
terials with the spatial dispersion it is more convenient 
to use the surface impedances instead of the dielectric 
functions . These impedances are connected with the 
dielectric functions by the relations 
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where the wave vector is k = (q^, Qy, kz) and the z-axis 
is perpendicular to the body surface. 

The rcfiection coefficients are expressed via the 
impedances as 

rs{q,uj) ^ -—f-j—^, rp{q,uj)^-— — (10) 



Lu/cko 



cko /uj + Zp' 



F{x) = 



dz- 



where ko = ^uP- j (? — q^ . Now we can use the refiection 
(6) coefficients pU)) in the Lifshitz formula. This formula is 
usually presented via the imaginary Matsubara frequen- 
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cies oj — !■ iQn = i2TTTn/H: 



T 



fl—S.p 71 — 
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where 



(11) 

= r^i{q,C,n) is the reflection cocfBcicnt of the 
body i {i — 1,2) for the polarization /i (/i = s,p). 

Let us consider the influence of the spatial dispersion 
on the force. When T is around room temperature typical 
values of the parameters are 7 ~ 10^^ — 10^^ rad/s, k ^ 
1/a, and v < 10^ m/s. Then the natural value of x in 
is large, \x\ ^ 1, for most of the materials in the 
interesting distance range 10-1000 nni. In this limit both 
dielectric functions become local: 



it = Si{k,iCn) = eo(«Cn) + 



Cn(Cn+7) 



n>0 (12) 



This is true, however, for n > 0. In the case n 
C ^ the result is different: 

,2 



et{k,iC) = £0 + 



(13) 



where £0 = £o(0). In the real frequency domain this limit 
is realized for <C 7. As one can see £; demonstrates the 
nonlocal character. 

We can conclude that at room temperature all the n > 
terms in ([TT|) can be calculated with the local dielectric 
functions. The n = term has to be corrected to take 
into account the nonlocality. Calculating the impedances 
([8|) and ^ with the functions (fT3|) and substituting them 
into the reflection coefficients (flUl) one finds 



0, 



£0 



(14) 



Here is zero because in the static limit s-polarized field 
is reduced to pure magnetic field, which penetrates the 
nonmagnetic material. One can see that Vp interpolates 
between a good metal and a pure dielectric. Really, im- 
portant values of q are q ~ 1/a, then for a <ti l]j = 
y/eo/ki) the refiection coefficient rp = {eq — l)/(£o + 1) 
is as for a dielectric with the permittivity Eq- In the 
opposite limit rp = 1 corresponds to good metals. 

In the local theory the n = term in the Lifshitz for- 
mula (llip can be presented in the Lifshitz form Q 
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where R = rpirp2 is just a constant. For example, for 
two dielectrics with permittivities £1 and £2 it is 



R 



(£1 - 1)(£2 - 1) 
(£l + l)(£2 + l)' 



(16) 



In the more general theory, which takes into account 
the spatial dispersion of the materials, the n = term 



can be presented in the same form (|15p but now R is 
not a constant but a function of q as follows from (|14p . 
Introducing for each material the parameter — 2a/ ho 
(i — 1, 2) we can present R = rpirp2 as 



R{x) 



■el 



(17) 



From IjlSP and (jl7p we can reproduce the main formula 
(34) in [13]. For that in the large distance limit we can 
consider the second body as rarefied £2 — 1 ^ 1, C2 — > 
and calculate the atom-body potential a,sV — — Fo(a, T). 

Our analysis of the Casimir-Lifshitz force for poor con- 
ductors at room temperature can be applied for the de- 
scription of a recent experiment [23l . |2J], the result of 
which did not found a reasonable explanation yet. In 
this experiment a p-type silicon membrane with the car- 



or rier density N ^ 5 ■ 10 



14 



cm 



was excited by the laser 



light. The density of photogenerated carriers was varied 
in the range Np^ = (1.4-2.1) • lO^^ cm.-^. The force dif- 
ference in the presence and in the absence of laser light 
was measured with an atomic force microscope (AFM) as 
a function of the distance between the membrane and the 
sphere at the end of the cantilever. It was found that the 
experimental results agree well with the dielectric mem- 
brane rather than with the semiconducting one. On this 
basis a controversial conclusion was made that one has 
to disregard the finite conductivity of silicon when the 
Casimir-Lifshitz force is evaluated. 

The controversy manifests itself in the local variant of 
the theory as discontinuity of the n = term in the Lif- 
shitz formula. For a dielectric (body 1) and metal (body 
2) this term is given by (|15p with R = {ei — l)/(£i + 1). 
Even for infinitely small conductivity of the dielectric this 
term has to be calculated with R — 1 and the result will 
coincide with the n — term for two metals. The nonlo- 
cality smooth out this discontinuity allowing continuous 
transition from dielectric to metal in accordance with the 
common sense. 

Figure [1] shows the experimental data for two densities 
of the photogenerated carriers Note that Nph in both 
cases are the same within the errors. Comparing figures 
[IJa) and[Hb) one can see that the experiment can hardly 
distinct between the nonlocal theory and the local theory 
with zero conductivity of Si. 

The controversies of the local theory are closely related 
with the behavior of the entropy at T ^ 0. Let us check 
the Nernst heat theorem in the nonlocal case. At low 
temperatures we cannot use the dielectric functions (jl2p 
and (fT3|) since the variable now is not large. Actually, 
when T is sufficiently small, the opposite limit is real- 
ized \x\ ^ 1, for which the nonlocal effects are strong. 
This is because ^{T) decreases with T faster than lin- 
early. Important imaginary frequencies contributing to 
the temperature dependent part of the free energy are 
C ~ 2'kT /h. Now the nonlocality is important for many 
terms in the sum (jlip . In this limit the dielectric func- 
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FIG. 1: (Color online) Difference of the forces on the low- 
doped Si membrane in the presence and in the absence of 
photogenerated carriers (dots with bars). The gray stripes 
are predictions of this work accounted for errors in Nph ■ The 
blue curves are predictions based on pure dielectric Si for the 
central values of Nph- (a) and (b) correspond to different laser 
powers exciting the carriers (see H^). 

tions are 

et(zC,fc) =eo(*C) + ^, ei{iC.k)^e,,{iO + ^- (18) 

Here we suppressed the index n and a is 37r/4 for the class 
I and -^7? for II. Note that 7 falls out from the result; its 
role plays the Landau damping frequency vk. 

Luckily there is no need to calculate the free energy 



with ((18)1 : all the work was done before. For the class I 
materials the calculations were performed in [T^ but the 
distance dependent leading term was presented in [20l |: 



T K -0.0193 



hvpc^ 



(19) 



The spatial dispersion changes the behavior of the free 
energy with T: instead of linear as in the local theory 
[3] it becomes quadratic. It ensures the right behavior of 
the entropy S ~ +T. 

Different situation is realized for the class II materials. 
In this case the density of carriers N{T) is a function 
of temperature. Due to presence of the energy gap for 
these materials at T — > the dependence is exponen- 
tial N{T) ~ e~^/"^, where A is the gap. Because both 
parameters uj^ and are proportional to N{T) their ef- 
fect in the dielectric functions is exponentially small. 
This conclusion is also true for the hopping mechanism of 
conductivity. In this case the effective density of charges, 
which are able to move, is ~ e~^^^ i25]. In the limit 
T ^ the class II materials become pure dielectric. For 
dielectrics the entropy disappears with T as S ^ T'^ [1] . 

In conclusion, the Lifshitz formula was generalized to 
materials, which cannot be described as local media. The 
spatial dispersion naturally resolves paradoxes appearing 
in the theory for local materials with finite conductivity. 

When this work was finished a preprint [26] appeared 
where the Lifshitz theory was generalized to semiconduc- 
tors. Ideologically the approach is similar to ours but less 
general and technically different. 

Fruitful discussion with L. P. Pitaevskii is appreciated. 
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